The Energy of a Trapped Interacting Bose Gas 
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A Bose gas in an external potential is studied by means of the semi-classical approximation. 
Analytical results are derived for the energy of an interacting Bose gas in a generic power-law 
trapping potential. An expression for the chemical potential below the critical temperature is also 
obtained. The theoretical results are in qualitative agreement with a recent energy measurement. 
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The condensation of an ideal Bose- Einstein gas is one of the most striking consequences of quantum statistics [Q. 
When Bose-Einstein condensation (BEC) occurs at a sufficiently low temperature, the zero momentum state can 
become macroscopically occupied. For many years, it was considered hopeless to experimentally observe BEC in an 
atomic gas with weak interactions. With the development of techniques to trap and cool atoms, BEC was recently 
observed directly in dilute atomic vapors |^-^ • The new experimental achievements have stimulated great interest in 
the theoretical study of inhomogeneous Bose gases. 

The thermodynamic properties of trapped atomic Bose gases undergoing BEC can be altered by the spatially varying 
trapping potential. The interaction between atoms may have a significant effect on the thermodynamic properties. 
There have been several investigations analyzing the dependence of the critical temperature on the trapping potential 
and weak interaction in the Bose gas |^-|^ . The thermodynamic properties of Bose gases in an external potential have 
also been discussed in Ref s. 

In a recent experiment ]10|], after turning off the trapping field, the kinetic energy of a sufficiently expanded atom 
cloud was measured. The experimental data is in a good agreement with the theory of a trapped ideal Bose gas 
when the temperature is above critical temperature. However, a discrepancy with the theory of ideal Bose gases is 
found below the critical temperature. To explain this discrepancy it is important to know how the mutual interaction 
affects the energy. Here we shall derive some analytical expressions for the energy of a trapped non-interacting and 
interacting Bose gas under the semi-classical approximation. 

First, let us calculate the kinetic and potential energies of an ideal Maxwell-Boltzmann gas and of an ideal Bose- 
Einstcin gas in an external potential. For an ideal Maxwell-Boltzmann gas, according to the Maxwell-Boltzmann 
distribution, the local atom number density is given by 



, (r, p, T) = exp[-/3(pV2m + V{y)- , 



(1) 



where (3 = kT and ^ is the chemical potential which comes from the normalization. Integration of the number density 
(|l|) over the whole phase space gives the total number of atoms 
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For the generic power-law potential discussed in Rcf. 

V{r) = £1 

where Li, L2 and L3 are the linear sizes of the volume, and £1, S2, £3, P, I and q the parameters, Eq. 
with 
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In the derivation we have used the formula for the gamma function 
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r(z)= / t'-^e-* dt. (7) 
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Relation (|^) determines the chemical potential /i from a given total number N of atoms. The total kinetic energy can 
be obtained as 



i^MB = / 1^ dN{r,p,T) = hr^ J exp[-/3y(r)] dr = ^NkT, 

which is just a result of the energy equipartition principle. Similarly, the total potential energy is 
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Thus, the total energy of the trapped ideal Maxwell-Boltzmann gas is 

El, = {^ + l)NkT. (10) 

We next consider an ideal Bose-Einstein gas in the generic power law potential. Under the semi-classcial approxi- 
mation |lT|J^j7[| the local number density is given by the Bose-Einstein distribution function 

neB(r, p, T) = {exp[/5(pV2m + ^(r) - /i)] - (11) 

from which the reduced spatial distribution of atoms is PJTl[] 

p(r) = A-3ff3/2 [exp(-/3(l/(r) - (12) 

where gu{x) = T,j=i I f ■ 

The critical teinperature for an ideal Bose-Einstein gas in the generic power law potential to undergo BEC has been 
found in Refs. §,§ to be 
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where (^{v) — 5i/(l) is the Riemann zeta function. 

The cases of the temperature below and above Tc should be treated separately. For T > Tc, integrating the number 
density function ( |ll] ) over the phase space, we obtain the total number of atoms j|] 

^ = /^"'/ nUr,P,T)dvdp^^^^g,+,{z). (14) 

In the derivation we have used the definition of the function f^^. Similarly, we have the total kinetic energy 

K^E = h-^J ^n^E{r,p,T)dpdr =^^^^0^ g,^+2iz), (15) 

and the total potential energy 

= / V{r)pir)dr = ^^M| + T + ^) ^^^^^'^^ (^6) 
The total energy of the trapped ideal Bose gas is then |l^ 

Ei^ = ^-^(v + l)NkT (17) 

Comparing with the trapped ideal Maxwell-Boltzmann gas, we see that the only difference in the two total energy 
expressions is the factor 17^+2 (2)/.'777+i(^) for the Bose gas, which breaks the linear dependence of the total energy on 
the temperature. 

When the temperature is below T^, the contribution to the energy of the system comes only from the normal 
component of the Bose gas. From Eqs. (p^, (|l5|), (16) and ( p7| ) the kinetic, potential and total energy of the system 
can be written as 
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The energy of the system will go to zero when the temperature approaches zero. 

So far we have not considered the role of the mutual interaction of Bose atoms. For a trapped interacting Bose 
gas, we divide the total energy into three parts: the kinetic energy Ki^t, potential energy Vint and interaction energy 
Uint, i.e. 



^int — ^int + Vint + Ui, 



(21) 



In the local density approximation or semi-classical approximation, the local number density function of a trapped 
interacting Bose gas is given by [^11 



i(r, p, T) = {exp[/3(pV2"i + V{r) + AaX'^p{r) - ^)] - 1}"^ 



(22) 



In similarity to the non-interacting Bose gas, the cases of the temperature above and below Tc need to be dealt 
with separately, where Tc is the critical temperature of the trapped interacting Bose gas S . At a temperature above 



Tc, from (|22|) the total number of atoms is 
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In the derivation we have expanded the exponential density function with respect to the small parameter a/ A, and 
kept only terms up to the lowest order in a/ A [||j9|. Up to this first order correction, Eq. ( p^ can be used. By the 
same procedure, the total kinetic and potential energy are obtained as 
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Similarly, with ^(r) being replaced by 2a\^ kT p{r) , the interaction energy is calculated as 
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From the above three expressions for Ki^t, Vint and Uint, we see that the correction to the energy of the trapped ideal 
Bose gas due to the mutual interaction of atoms is proportional to a/ A, hence is very small for a weak interaction. 
This explains why the ideal Bose gas result agrees well with the experimental data above Tc. 

When the temperature is below Tc, a region where a great number of atoms are in the BEC state forms. In the 
local density approximation, a relation between density and potential inside the region reads |0] 



V{r) + 47rap(r)?i^/m = V{ro) + inapoh^/m, 



r < ro. 



(28) 



where ro corresponds to the edge of the region. At r = Tq", we have pir^) = po — A^^(73/2(l), which is the critical 
density of the gaseous phase for the untrapped homogeneous ideal Bose gas at temperature T. From the above 
relation, we find the density of the condensed atoms as the difference of p{r) and po 



Ps{r) = p{r) - po = 
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In order to obtain ps, the potential V{ro) at the edge of the condensed region needs to be determined. Integrating ps 
over the whole volume, we have the total number of atoms in the condensed state 



mx I{p,l,q) . 



where the derivation involves the integral 
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At the weak interaction limit we may approximate iV^ and Tc by those for the trapped ideal Bose gas [plql, i.e. 
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where we have written Tc as Tc- From Eqs. (^) and (|3^) we find 
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For the condensation region at T < there exists another relation [ pA| 
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Setting r = Tq in this formula, from Eq. ( p3|) wc find the chemical potential as a descending function of temperature 
to be 
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which, for a cylindrically symmetric harmonic trapping potential, reduces to 
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where uq is the geometric mean characteristic length of the trap, i.e. Oq = h/ (mwo) with luq = uj\uJz. This is consistent 
with the Ginzburg-Pitaevskii-Gross mean field theory |]l3| . 

According to Refs. ||ll|,|l^, the density of the mutual interaction energy in the condensation region is 



u^nt = 2a\^kTp^{v) - a\^kTpl{v) 

= a\^kT[2pl{v) + Apopsiv) + pI{v)], (r < ro). 



which, after integrating over the condensation region (r < rp), gives 
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The first two terms depending on po vanish when the temperature approaches zero. So, at low temperature the main 
contribution to Uint comes from the third term corresponding to the interaction within the condensed atoms. We 
may then neglect the first and second terms of (|3^) to write 
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To an approximation, we may neglect any correction to the energy from the normal component, and express the 
energy of the trapped interacting Bose gas as that of the trapped ideal Bose gas plus the only correction from the 
mutual interaction within the condensed component. 
In the experiment the trapping potential is [|lO| 



V{v) = -muiri + -uiri, (40) 

which is harmonic. In the ideal Bose gas approximation, at a temperature T above T^, Eqs. (^5|) and ( |l6| ) reduce to 

334(2:) 
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where the fugacity z is related to the total number of atoms as 
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which yields 53(2) = {Tc/TY'gi{\). The energy of the atom interaction is negligible at T > Tc- 

When the temperature is below Tc, from Eqs. (Iq) and (p^), the approximate kinetic and potential energy are 



2C(3) \T, 
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(43) 



The contribution to the energy from the interaction within the condensation component is now significant. From 
Eqs. ([l3| ) and (|39|), this interaction energy is 
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In the experiment the trapped atoms are ^^Rb ||T^. The s-wave scattering length of ^^Rb is a « 100 in the unit of 
the Bohr radius as ■ The trapping harmonic potential is axially symmetric with the frequency ratio of axial to radial 
being VS- The final stage of the evaporative cooling is performed at i/^ = 373 Hz. The trapping field is then turned 
off to initiate the expansion of the atom cloud non-adiabatically. The kinetic energy iJ^^p of atoms is measured after 
a sufficient expansion which ensures the interaction energy to convert to purely kinetic energy. Therefore, we may 
estimate E„^^ as the sum of and , i.e. 
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for T < Tc, 



(45) 
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Assume that the total number of atoms in the trap is 40000. The scaled energy per particle, E/NkTc, is plotted 
versus the scaled temperature, T/Tc, in Fig. 1, where the solid, dashed and dotted lines correspond to the ideal Bose 
gas, ideal Maxwell-Boltzmann and interacting Bose gas, respectively. When the temperature is rather high, they are 
very close. However, at low temperature the differences among them are prominent. The theoretical result of the 
interacting Bose gas is in qualitative agreement with the experimental data. We have considered only the contribution 
from the mutual interaction of the condensed atoms. If the finite size effect |l5| , p^ and all other interaction terms are 
included, a quantitative comparison with experimental data can be made. 
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Figure Captions 

Fig. 1 Plot of the scaled energy per particle E/NkT^ of trapped gas vs. scaled temperature T/T^. The solid, dashed 
and dotted lines correspond to an ideal Bose gas, ideal Maxwell-Boltzmann gas and interacting Bose gas, respectively. 
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